Exact Results on Dynamical Decoupling by n Pulses in Quantum Information 

Processes 
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The aim of dynamical decoupling consists in the suppression of decoherence by appropriate co- 
herent control of a quantum register. Effectively, the interaction with the environment is reduced. 
In particular, a sequence of n pulses is considered. Here we present exact results on the suppression 
of the coupling of a quantum bit to its environment by optimized sequences of n pulses. The effect 
of various cutoffs of the spectral density of the environment is investigated. As a result we show 
that the harder the cutoff is the better an optimized pulse sequence can deal with it. For cutoffs 
which are neither completely hard nor very soft we advocate iterated optimized sequences. 

PACS numbers: 03.67.Pp,03.67.Lx,03.65.Yz,03.65.Vf 



I. INTRODUCTION 

Almost six decades ago in 1950 Hahn demonstrated 
that spin echos in liquid NMR can be obtained by apply- 
ing a 7r pulse in the middle of a time interval [l|. This 
idea was developed further by Carr and Purcell who pro- 
posed iterated cycles of two ir pulses to reduce the effect 
of unwanted interactions Q. Further refinements were 
introduced by Meiboom and Gill [|[ . Since then this tech- 
nique of coherent control has been established in NMR, 
see e.g. Ref. 0. 

The fascinating possibilities of quantum information 
have stimulated a great interest in the coherent control 
of small quantum systems, see e.g. Ref. i. The idea to 
preserve coherence by iterated 7r pulses periodic in time 
was rediscovered in the context of quantum information 
by Viola and Lloyd [6| and by Ban [7| in 1998 for a spin- 
boson model and subsequently generalized to open sys- 
tems [§[ ; a short review is found in Ref. H|. For symmetry 
groups with inefficient representations randomized pro- 
tocols are advocated, see Ref. [1(3 and references therein. 

Recently, periodically iterated Carr-Purcell cycles have 
been advocated for the preservation of the coherence of 
the electronic spin in quantum dots [Til . Il2j . Besides 
periodic iteration of pulse cycles also concatenations of 
cycles were proposed and it was shown that they suppress 
decoherence in hi gher orders t in the length of the time 
interval t [ll|, fl3l Il4, [H[ • But the achieved exponent I 
grows only logarithmically in the number of pulses n. 

In parallel, the author showed that neither the itera- 
tion nor the concatenation of the Carr-Purcell two-pulse 
cycle is the optimum strategy for a single-axis bosonic 
bath. [l6j. Cycles with n pulses at the instants Sjt 



[nj/{2n + 2)] 



(1) 



achieve the optimum suppression of decoherence in the 
sense that any deviation of the signal occurs with a high 
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power in t, namely t 2n+2 where n is the number of pulses, 
i.e. only a linear effort is required. The Carr-Purcell cycle 
is retrieved for n = 2. Up to n < 5 the result ([I]) was 
previously shown for general models [13] • 

Then, Lee et al. [IH observed in numerical simulations 
of spin baths that the pulse sequences obeying Eq. ((T|) 
suppress also the decoherence for spin baths. They could 
also show analytically that the sequence defined by (Q]) 
works for the most general phase decoherence model up 
to n = 9. We have been able to extend this analytical 
proof up to n < 14. An unrestricted derivation, however, 
is still lacking. Lee et al. [HI] also argued that the opti- 
mized sequence ([l} works well only when the expansion 
in time is applicable. 

There is a multitude of experimental results in NMR on 
coherent control and the suppression of decoherence, see 
e.g. Ref. |4|. We highlight results in the context of quan- 
tum information related to pulse sequences [l9| : for a 
overview see Ref. |2fj. But also in semi-conductor physics 
there are many encouraging results in pro l onging the co- 
herence time of a qubit by n pulses [2ll.l22l 1 2 31124, l2o| . In 
experiment, one must trade off between the advantages of 
the suppression of decoherence by multiply applied pulses 
with the detrimental effects of imperfect realizations of 
pulses, for instance the finite duration of a pulse so that 
it cannot be regarded as instantaneous [2t3 ]. 

The aim of the present article is threefold. First, we 
provide the explicit calculations leading to the important 
relation (fT]). Second, we generalize the previous result 
[il| on a particular signal to a statement on the unitary 
time evolution. Thereby, we provide the general proof 
for the applicability of |T]) for an arbitrary initial quan- 
tum state. Third, we use various spectral densities J{ui) 
in the spin-boson model to discuss under which condi- 
tions the optimized sequence works well, namely when 
the high-energy cutoff of the decohering environment is 
hard enough. To cope with medium hard cutoffs we 
propose iterated sequences of short optimized cycles of 
pulses. 

The article's setup is as follows. In the following Sect. 
|TT]the explicit calculation for the spin-boson model is pre- 
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sented, both for the signal in a generic decoherence ex- 
pcrimenet and for the general time evolution. The results 
are also given for classical noise. The subsequent Sect. Mil 
treats the general phase decoherence model. Section IIVI 
presents a discussion of the applicability of the optimized 
sequences and establishes a link to the nature of the high 
energy cutoff. The conclusions IVl summarize the results 
and discuss their implications for further developments. 



II. SPIN-BOSON MODEL 

We consider the model given by the Hamilton operator 
H = Y J ^b\b l + l -a z Y J Hb\ + b l ) + E (2) 

i i 

consisting of a single qubit interpreted as spin S = 1/2, 
whose operators are the Pauli matrices a x ,a y , and a z . 
The environment is given by the bosonic bath with an- 
nihilation (creation) operators . The constant E sets 
the energy offset. The properties of the bath are defined 
by the set of real parameters {Ai, uii\. This information 
is conveniently encoded in the spectral density [13, HH 



J(u)=J2\ 2 i S(cj-u> i ). 



(3) 



Obviously, H in ([2]) does not allow for spin flips since it 
commutes with a z . Physically this means that the decay 
time T\ of a magnetization along z is infinite. But the 
decoherence of a magnetization in the cr^cTy-plane is cap- 
tured by H so that the decay time Ti can be investigated 
in the framework of this model. 

The Hamiltonian H in @ is analytically diagonaliz- 
able. For any operator A we will use the notation 



A eS := UAUK 



(4) 



The unitary transformation U is chosen so that it diago- 
nalizes H 



(5) 



The appropriate unitary transformation is 

U = exp(a z K). (6) 
The operator K is anti-Hcrmitean 



so that U is indeed unitary. The energy offset after the 
transformation reads 



AE = E- 



duj. 



(8) 



But the global energy offset is not measurable, so that 
its quantitative form does not matter. 



A. Signal without n Pulses 

Here we discuss the simple experimental setup without 
any 7r pulses. We start from the state | "[")■ Then a 7r/2 
pulse is applied to rotate the spin from the ^-direction to 
the a;y-plane. To be specific, we rotate the spin about x 
by the angle 7 with the help of the unitary transformation 

As (7) : = exp(-i-ya x /2) (9a) 
= cos(7/2)+icr a! sia(7/2). (9b) 

The rotation is best seen by stating that 

D x (j)^<7 z D x ('j) = a z cos 7 + a y wa.j. (10) 

For 7 = 7r/2 a spin along z is turned into a spin along y. 
We will use D x (n/2) = (1 + ia x )/y/2. 

In the xy plane the spin will evolve. After the time t 
a measurement of <r y yields the signal 

s(t) = (t |A c (vr/2) t cxp(^i^^) f 7 ^/ exp(-^i^^)A c (7r/2)| t). 

(11) 

Since H does not induce spin flips and (f \a y \ f) = = {[ 
\(T y \ I) we know that (f | exp(iHt)a y exp(— iHt)\ f) = 
and (| \cr x exp(iHt)cr y exp(— iHt)cr x \ f) = 0. Hence the 
signal is given by 

s(t) = 3m(t \a x exp(iHt)(T y exp(-iHt)\ |). (12) 

To evaluate this expression explicitly we change to the 
basis in which H is diagonal 



s(t) = 3m(T \af eMiH ea t)a e y n eM-iH t)\ T>- (13) 

Note that the state | f) is not altered by U. For the 
explicit calculation of the effective operators we use 



-over 



x/y 



(14) 



and obtain 



<jf/ y = exp((T z K )a x / y exp(-<7 z K) (15a) 
= exp(2(T z K)a x / y . (15b) 

Hence the action on particular spin states is 

af IT/I) = exp(T2A-)U/t) (16a) 
af\]/[) = ±zexp( T 2if)| 4/T), (16b) 

where either the first spin orientation and the upper sign 
holds or the second spin orientation and the lower sign. 

Turning to the time dependence we define generally 
the time dependent operators 



A(t) := exp(iH cS t)Aeiq>{-iH eS t). 



(17) 



Note that H contains only the bosonic degrees of free- 
dom and it is diagonal. Hence it is easy to see that 



b\(t) = &}exp(tw<i) 
= h exp(-iwit) 



(18a) 
(18b) 
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whence 



whence we obtain 



K(t) = y exp{iu>it) - b t expC-iwit)). (19) 



x( 



(i)=Eril 4sin2 (^*/ 2 )^+ 6 ^l) 



(30) 



With these definitions the identities (|T6|) apply also to 
the time dependent operators cr cfI (t) and K(t) 



<C(*)IT/i) 
<(*)IT/i> 



exp(T2Jf(t))| 4 / T> (20a) 
., .. , , - ±iexp(T2A-(t))U/T>. (20b) 

With these identities we can write for the signal 

s(t) = 3m(T \af(0)af(t)\ T) (21a) 
= 3mi(l |cxp(2if(0))exp(-2if(i))| |) (21b) 
= !Ke(exp(2i<r(0))exp(-2.fr (*))), (21c) 

where we took the expectation value in the spin sector 
in (|21c|) so that only a bosonic expectation value with 
respect to the bilinear H cS must be computed. This is 
eased by the Baker-Campbell-Hausdorff (BCH) formula 

exp(A) exp(B) = exp(A + B) exp(L4, B]/2) (22) 

which is valid if [A, B] commutes with A and B. This 
yields 

s(t) = ^eexp(-2[^(0),X(t)])(exp(-2Ais:)> (23) 

with AK := K(t) — K(0). Any expectation value of an 
exponential of a linear bosonic operator A with respect 
to a bilinear Hamiltonian such as H can be reduced to 
the exponential of an expectation value by 



(exp(A)) = exp((^ 2 )/2). 

Hence we have 

s(t) = meexp(-2[K(0),K(t)])exp(2(AK 2 )) 
which simplifies due to the Hermitecity of AK 2 to 

s(t) = cos(2<p(i)) exp(-2x(t)) 
where the phase is given by 

<p{t) := i[K(0),K(t)} 



A? 



2uj 2 K ' 



— — sux(ojt)au>. 



The exponential suppression is given by 

'2\ 



x(t) 



-(AK- 



where 



AK 



•iOJit 



(24) 

(25) 

(26) 

(27a) 
(27b) 

(27c) 

(27d) 

(28) 
1)1 (29) 



The bosonic occupation is such that the last expectation 
value equals coth(/3wi/2) so that we finally have 



x(t) 



sin 2 (wt/2) 



coth(/3w/2)dw. (31) 



This concludes the derivation of the signal without any 
dynamical decoupling. The formulae (6) in Ref. [l6| are 
rederived in all detail. The above derivation sets the 
stage for the derivation in the case of dynamical decou- 
pling by sequences of 7r pulses. 



B. Signal with n Pulses 

Here we consider a sequence of tt pulses which are ap- 
plied at the instants of time Sit with i G {1, 2, . . . , n} so 
that n pulses are applied and the total time interval t 
is divided into n + 1 subintervals. For notational conve- 
nience we set Sq = and <5 n +i = 1. It is understood that 
5 i+ i > Si for all i G {0, 1, 2, ... , n}. 

The 7r pulses are taken to be ideal, that means they 
are instantaneous so that during their application no cou- 
pling to the bath needs to be considered. The possible 
workarounds if this is not justified in experiment are dis- 
cussed elsewhere [26[ . For simplicity, we take the ir pulses 
to be realized as rotations about a y 

D y (j) := exp(-ija y /2) (32a) 
= cos(7/2)+io-„sin(7/2). (32b) 

which implies for 7 = n simply D v (tt) — ia y . Below we 
will use a y only because the factor i corresponds to an 
irrelevant global phase shift. 

The signal is given in general as before by 



»(*) = 3m(T \<T x tf<TyR\ t) 



(33) 



where the time evolution is changed its form exp(— iHt) 
in |H to 



R 



e -mSn + ,-S n )t (7ye -iH(S n -S n ^)t (7v 
e -iH(S3-Si)t -iH(Si-S )t 



(34) 



The expression (|34[) becomes much more compact if writ- 
ten after the diagonalisation by U as given in © 



R cti = e' lH *i? off 



(35) 



where we express R based on (fTT)) in the form 

R cS = ofVnt) af(S n ^t) . . . af(S 2 t) of (36) 
Then we arrive easily at 

s(t) = Jm(T \af (0)R° S ^af(t)R eS \ T) . (37) 
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Equation (|37p can be converted by means of (|20|) into 
the following purely bosonic expression 

s(t) = 3m( e 2 ^ > ie~ 2K ^ (-i)e 2K ^ . . . 

... (-i)(-l)» e C-l) B 3*<*»*) 

(-i)(-l)«+ 1 e (- 1 ) n+l2 ^W 

(-ix-l)^- 1 )" 2 ^*-*) ... 

...(-t)e 2Jf ( &i )i e - w < w A, (38) 

Counting the factors (—1) and i one finds that they all 
combine to a single factor i. This is easiest seen by 
combining the prefactors in front of each of the terms 
(— iy i e^~ l ' 2Kl ^ Snt ' which all occur twice so that each 
kind of these terms provides a factor (— i) 2 = — 1 yielding 
a total factor (-1)™. This is multiplied with (-i)(-l) n+1 
from the prefactor of e^ -1 )" 2 ^ 15 "*' which is the only term 
occuring only once. Hence we have 



S(t) = me ^K(0) e -2K(S lt ) e 2K(S 2 t) 



e (-l)"2K(S„t) e (-ir +1 2if(t) e (-l)"2K(5„t) 
_ e 2K(S 2 t) e -2if(5!t) 



(39) 



Applying the BCH formula ([H]) yields 



s(t) = m.texp{2i<p n (t)) (exp(2A n K)) (40a) 
= cos(2p n (i))exp(-2xn(t)) (40b) 



where we used the identity (|24p to obtain the second line. 
Therein the suppression Xn(t) '■= — (A n K 2 ) results from 



A n K := K(0) + (-l) n+1 K(t) + 2Y,(-l) n K(S t t) (41a) 

i=l 

= T,£Mv<>fa t )- b iVnW)> ( 41b ) 



where 



l {z) :=l + (-l)" + V z + 2^(-l)V^. (42) 



Thereby we arrive at 



W = £ ^2 l2/n(^*)| 2 (&& + Ml) (43a) 

j( w )lM4^coth(/3w/2)dw. (43b) 
o ^ 

The phase (f n (t) in (|40p can easily be computed by 
the following trick. Using (|22[) we combine the sec- 
ond and third factor in (|39p . i.e., the two exponen- 
tials e ~ 2K ( s i t ) e 2K(S 2 t) ^ ^ Q Qne ex p 0nen tial and the last 
and last-but-one factor, i.e., e 2K(52t) e - 2K(Slt K Obvi- 
ously, the occuring phases cancel. This procedure can 



be repeated by including the factor e~ 2K ^ 3t ^ next both 
in the growing last exponential and in the second ex- 
ponential. Iteration up to and including the factor 
e (-i) 2K(s n+1 t)^ winch can be thought as being split 

into e (-ir + 1 2if(<Wl*)/2 e (-l)" + 1 2if( 5 „ +lt )/2 ; leadg tQ twQ 

exponentials whose respective arguments contain all term 
K(6jt) except the very first K(0). Furthermore, the two 
respective arguments are equal so that the exponentials 
can be combined without further phase yielding 



s(t)=mc(e 2K ^e 2A " K ~ 2K (°A 



(44) 



From this equation we arrive at (|40a[) by defining 

<p n (t) := -i[K(0),A n K] (45a) 



J(w) 



(y n {iut) - y* n (Lut))dtu (45c) 
(45d) 



4w 2 

-—r-X n {Ujt)duJ 

2lo a 



where we used (|41bj) in the second line (|45b|) . The last 
line (145d|) reproduces Eq. (8b) in Ref. [l6| with 



x n (z) := i(y n {z) - y*( z ))/ 2 
= -%„(z) 



(46a) 
(46b) 



= (-l)"sin(z) + 2^(-iy +1 sm(z^), (46c) 

3=1 

where the last line (|46cp corrects Eq. (9) in Ref. [l6| in the 
factor 2 in front of the sum. 

Thereby, we have derived all the results used in the 
analysis in the previous paper [l6l |. 



C. Optimization of the Sequence of n Pulses 



A particular asset of the equations (|40bl43l45|) to- 
gether with (|46a[) is that it is obvious that any devia- 
tion of the signal s(t) from unity is kept as low as pos- 
sible if |yn(2)| is kept as small as possible. Note that 
this strategy holds equally well for (p n and for Xn- If 
y n is of the order p in some small parameter p, for in- 
stance p = t n+1 , then Xn = C(p 2 ) whence we deduce 
cxp(— 2xn) = 1 — 0(p 2 ). In analogy, we find cp n = 0{p) 
whence we deduce cos(— 2<p n ) = 1 — 0{p 2 ) so that both 
factors are close to unity in the same way. Hence the total 
signal s(t) is close to unity in this order s(t) = 1 — 0(p 2 ). 

So our aim is to choose the n instants {Sj} such that 
y n {z) is as small as possible. The best way to do so is to 
make the first n derivatives of y n (z) vanish. Note that 
2/n(0) = for any sequence {Sj}. The mth derivative 
reads (m > 0) 



d'; 



\n+l 



(47) 
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Hence we have to solve the set of nonlinear equations 

n 

= (-l) n+1 +2J2(- I ) j $7 ( 48 ) 
i=i 

for m € {1, 2, . . . , nj. For finite n, solutions can easily be 
found analytically [17| and numerically. Closer inspec- 
tion of these numerical solutions reveals that they are 
excellently described by the condition ([1]). 

Indeed, we can prove that ([!} is a valid solution for the 
set of equations (|4"5)) . To do so we choose a little detour 
by considering 



y™( h )\h= 



z/2 



exp(-iz/2)y n (z) 



(49) 



Obviously the equivalence 

y n (z) = 0(z n+1 ) & y n (h) = 0(h n+1 ) (50) 

holds so that the vanishing of the first n derivatives 
of y n (z) is equivalent to the vanishing of the first n 
derivatives of y n {h). The choice ([T]) implies by standard 
trigonometric identities 



5j = l/2-cos(7rj/(n + l))/2. 
Inserting this choice into y n (t) yields 



y n (h) = e- th + (~l) n+1 e lh + 



(51) 



-2^(-l) J cxp[-i/icos(7rj/(n+ 1))] (52a) 



= (-l) J exp[-i/icos(7rj/(n+ 1))]. (52b) 

j=—n—l 

Obviously y n (0) = 0. The mth derivative (m > 0) reads 



since the denominator in (|55ap remains finite in this 
range. Hence d m — and we know j/„(/j) = 0(h n ) and 
hence y n {z) = 0{z n ). This concludes the formal proof 
that fl]) represents a valid solution of the set of nonlin- 
ear equations (|48|) . We have not presented a proof that 
this is the only solution. But we presume that it is the 
only one which is physically meaningful with consecutive 
values 5j + i > Sj. 



D. Classical Noise with ty Pulses 

In Ref. [TH we argued that the fact that the optimized 
sequence (TTJ) works independently from the precise tem- 
perature indicates that it applies also to classical, Gaus- 
sian noise. The argument runs qualitatively as follows. 
Because ([TJ is the optimum sequence for all temperatures 
it holds of course also for T — ► oo. In this limit, the ther- 
mal fluctuations dominate over all the quantum effects 
and the bath behaves completely classically. 

A crucial corollary is that the pulse sequence can be 
used for all kinds of bath at elevated temperatures be- 
cause all physical systems behave like classical, Gaussian 
baths at high temperatures. Hence the applicability ex- 
tends beyond the spin-boson model discussed so far. Wc 
will discuss the general validity of (fTJ) in more detail in 
the next section. 

Here we present the calculation for classical noise in 
order to establish a quantitative relation. We consider 
the decoherence due to 



H = f(t)a z 



(56) 



where fit) is a random variable with Gaussian distribu- 
tion. [32| It is characterized by the expectation values 



d hVn\ h= , = {^) m £ (-l)''a* ro (7rj/(n+l)). (53) 

j=—n—l 



We compute explicitly d m :— (2i) m d™y n \ h _ 

n 

d m = £^ ^\y^ e ™j/{n+l) _ e ~inj/(n+l)yn 
j=—n—l 

m \ ( i-Kj{2v — m + 1) 



exp 



n+1 



V— j= — 71 — 1 

The last sum, however, vanishes for m < n + 1 
/ mj{2u — to + 1) N 



ex P 



3=-n— 1 



77 + 1 



(54a) 
.(54b) 



/ -, r +i exp(-77r(2;/ - m)) - exp(z7r(2^ - 777)) 
1 ' l + exp((*7r(2i/-m))/(n + l)) 1 ' 

= (55b) 



(/(*)) - (57a) 
(/(ti)/(*a)> = 9(h-h). (57b) 

Note the translational invariance in time. Then the signal 
s(t) after a tt/2 pulse reads 

s(t) = (T |£» x (7r/2)t e iF W- a y e~ iF ^ D x (n/2)\ t) 
= Jm(t \a x e iF ®' T *cT y e- iF W' T '\ T) (58) 

where F(t) := J Q f(t')dt' is the primitive of f(t). Since 
o x .y only flip the spin, see for instance Eq. (fTS)) for K = 0, 
we may write 



s(t) = ( e -^(*)e- iF «) 

= e -2(F(t) 2 ) 



(59a) 
(59b) 



where we exploited the properties of Gaussian random 
variables to obtain the second line (|59b[) . The exponent 
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can be computed easily 



(F(t) 2 ) = f dt t f dt 2 (f(h)f(t 2 )} (60a) 
Jo Jo 

(60b) 



= 2 [ dtx [ dt'g{t') 
Jo Jo 



I P[ ^lsm 2 (ujt/2)duj. (60c) 

1* 



7T In W 



g(t) = — p{u) cos(ujt)duj 



where we used g(t') — g(—tf) to obtain (|60b[) and the 
Fourier representation for (|60cj) 



(61) 



based on the power spectrum p(w). The comparison with 
the quantum mechanical result (|26I31[) yields exactly the 
same form except that tp(t) = because there are no op- 
erators which might not commute with themselves. The 
argument of the exponential %(t) = (F(t) 2 ) is identical if 
we identify 



(4/tt)pH = J(w)coth(/3w/2). 



(62) 



This provides the quantitative correspondence between 
the classical calculation and the general quantum me- 
chanical one. 

The extension to the signal in presence of the n pulses 
is also straightforward. The signal is given as in (155)) 
except that the time evolution R c \ is classically given by 



e -*- f W dt a 



(//, 



(63) 



Again, the dynamics of the spin is easily computed since 
it flips at each a y or a x according to (fTtj|) . The final result 

is s(t) — e~ 2 ( Fn< & ) where F n (t) is given by 



F n (t) := / f(t') Sn (t')dt' 
where s n (t l ) switches the sign according to 



(64) 



for t' < 
s n (t') := { (-1)3 for Sjt < f < S j+1 t (65) 
for t' > t 

for j € {0, 1,2,..., n}. Note that the Fourier transform 
s„(lu) of s n (t') is given essentially by y n (u)t) 



s n (t') exp(iujt')dt' = — y n (ut). (66) 

oo W 

Next (F n (t) 2 ) is expressed as convolution and integral 

/•oo />oo 

(F n (t) 2 ) = / / dt 1 dt 2 s n (t 1 )g(t 1 -t 2 )s n (t2) (67a) 



oo J —oo 



1 



p{w) 



\y n (ut)\ — ~-du. 



(67b) 



For the last line, Fourier transformation, Parseval iden- 
tity and the symmetry of the integrand are used. Again, 
we retrieve the quantum mechanical result (|40b|43bj) ex- 
cept for the phase ip n (t) which does not occur at all in 
the classical framework. The necessary identification is 
the same as before ([621) . 

We conclude that the classical decoherence and the 
one due to a quantum bosonic bath coincide except for 
the phases if the power spectrum 4p(w)/7r is identified 
with the product of spectral density J(u>) and bosonic 
occupation factor coth(/3o>/2). Hence the optimization 
of the quantum model applies equally to the classical 
problem. Therefore, the optimization JT]) applies to all 
models with (commuting) Gaussian fluctuations. 



E. Unitary Time Evolution with 7r Pulses 

So far we focused on the signal sit) as it results from 
a measurement of a y after a ir/2 pulse around a x . This 
appears to be a special choice. But in view of the spin 
rotational symmetry about the z axis it is sufficiently 
general to guarantee that the coherence of an arbitrary 
initial state is preserved by the optimized pulse sequence. 
To corroborate this point and to prepare for the discus- 
sion of the most general model for phase coherence we 
discuss the time evolution operator R of the spin-boson 
model in this section. 

The unitary operator R is defined in (13^1) . Using the 
identities (|4I35[) we get 



R = Uh- iH ° l R ca U 



cHV 



R 



e- lH 4 R with 



(68a) 
(68b) 
(68c) 



Inserting the time dependent version of (I15b|) (cf. also 
Eq. (I20bl ) 



af(t) = exp(2a z K(t)) a y 
into (|3"6")) and using (|14[) yields for n even 

^1 = R -o z K{t) e 2cr M K{S n t) e ~2a z K(8 n ^ 1 t) 

I n oven 

e 2<TzK(5 2 t) e -2a,K(S!t) e <7 z K(0) 

while for n odd we arrive at 
R 



n odd y 



_ e 2a z K(S 2 t) e -2<r z K(6it) e <r z K(0) 



(69) 



(70) 



(71) 



These results are combined to yield for the total time 
evolution 



R 



1 



exp(-iH cS t) exp(-i<p n (t)) exp(a z A n K) 

(72) 

where the upper entry in the curly brackets refers to n 
even, the lower one to n odd. The multiple difference is 
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defined and computed in (f4Tj) . Combining all the expo- 
nents to a single one makes a phase 4> n (t) occur which 
can be computed by commuting the various expression 
K(8jt) as required in (|22[) . We do not give the explicit 
expression because we do not need it here. What is im- 
portant is that this phase is a global one. It is just a real 
number and it does not depend on the spin; no Pauli 
matrix occurs because o\ = 1. Similarly, H does not 
depend on the spin. 

To assess to which extent the time evolution depends 
on the spin state we consider the difference between the 
evolution of an f and of a j state. We define for n even 



R\ 



(T \R\ T) S pin 
(I \R\ Ospin 



and for n odd 



i? T := \R\ T) sp i„ 

Ri ■= *<T \R\ i)s P i„ 



(73a) 
(73b) 

(74a) 
(74b) 



where the subscript sp ; n signifies that we compute the 
expectation value only with respect to the Hilbert space 
of the spin. The bosonic operators remain unaltered. 
Then we consider 



A(t) := i? T -i?x 



iH ° { h e -i Vn (t) [gA„K 



(75a) 
- A " K ] (75b) 



as proposed by Lee et al. [l8| . From the last formula and 
(|41bp it is obvious that the influence of the spin state is 
small for general sets {cjj, A^} if and only if y n {z) is small. 
Quantitatively, one has 



y n (z)=0(z n+l ) «■ A(t) = 0(t n+L ) 



(76) 



Thereby, we have shown explicitly that the condition 
y n (z) — 0(z n+1 ) implies generally that the coupling be- 
tween any spin state, i.e., any state of the quantum bit, 
and the bosonic bath is efficiently suppressed if the pulse 
sequence obeys (HJ. Note that this holds for all choices 
of {uJi,Xi}. 



III. GENERAL QUANTUM BATH 

So far we considered the spin-boson model ([2]). One 
might think that the optimized sequence (Q} is useful 
only for this model [lij ]. This is not the case. 

The first evidence for the general applicability of ((T|) 
is the fact that classical Gaussian noise can equally well 
be suppressed, see Subsect. Ill Dt Conventional wisdom 
has it that any generic model with fluctuations will dis- 
play Gaussian fluctuations in its high temperature limit. 
If this is true the optimized sequence Jl} is applicable 
generally for high temperatures. Note that the "high" 
temperatures need not be really high. The inter-spin 
coupling of nuclear spins is so low that already 1 Kelvin 



is sufficient to put a system of nuclear spins at high tem- 
peratures. 

The second evidence was found by Lee et al. [l8| . They 
observed analytically for up to n = 9 that an expansion 
of A(t) in powers of t for a general model yields vanishing 
coefficients for the optimized sequence (TT]). On the basis 
of this observation they conjectured that the optimized 
sequence (fTJ) is generally applicable for the generic model 
for phase decoherence, also called single axis decoherence 
model 



H = a z A 1 +A , 



(77) 



where Aq and A\ contain only operators from the bath. 
Below we use the notation H± = ±a z Ai + Aq. 

This model does not include spin flips; hence it implies 
an infinite life time T\ as a completely general decoher- 
ence model would do. But the phase decoherence of a 
precessing spin in the xy plane is described in full gen- 
erality because we do not specify for which operator A\ 
stands and the bath dynamics is fully unspecified. It is 
described by Aq. Such a model is experimentally very 
well justified as the effective model in the limit of a large 
applied magnetic field which implies that other couplings 
between the quantum bit spin and the bath are averaged 
out, see for instance Refs. [29| andl30l. 

We investigate the time evolution R from to t with 
tt pulses at the instants 5jt where j S {1,2, ... ,n}. The 
tt pulses are assumed to be ideal; they are given by a y so 
that R is given again by (fM]) . Next, using (Til)) , we shift 
all the factors a y to the very left side yielding 



R _ J 1 I e -iH(_i)» (<S„ + i-<5 n )t e -iH ( _ 1)n _i (i5 n — (5„_i)t 



e -iH_(« 2 -«0i e -iH+(Si-S )t 



(78) 



where the upper entry between curly brackets refers to 
an even number n of pulses and the lower one to an odd 
number. 

We define the unitary operators U p as the product of 
the p + 1 rightmost factors on the right side of Eq. (|78|) . 
that means for < p < n 



U p (t) ■— e ~ i - ff (-i) p ( 5 p+ 1 ~ 5 p) t g^^t-i)?- 1 (,S p -S p -i)t 



e -iff_(«a-4i)t e -iH+(Si-S )t 



(79) 



This operator can be expanded in a Taylor expansion 
with coefficients C^r 

oc 

Upit) = ^ ^ &]r'C^rA mj A mj _ x . . .A m . 2 A mi . 

j=0 m£Bj 

(80) 

The set Bj contains all binary words m with j letters, i.e., 
rrii G {0,1} where rrii is the ith letter, i S {1,2, . . . , j}. 
Note that also leading zeros count. We use \m\ for the 
checksum of m, i.e., the sum over all letters |m| := 
2i=i m i- The number of letters j of m shall be denoted 
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by ||z?l||- Using B as the union of all Bj with j > we 
may denote the expansion by 



u P (t) = V (-it) 



meB 



\C—A 

Sir, 



A m2 A mi . (81) 



Obviously, the coefficients which matter in the end are 
those for p — n. The statement A(t) = 0(t n+1 ) corre- 
sponding to (|76p for the spin-boson model is equivalent 
to the vanishing of all the coeffients which are prefactors 
of terms depending on the spin state. This means that 
all Cfr with |m| odd have to vanish as long as n > | \m\ \ . 

So far no general proof is available that these coeffi- 
cients vanish for the sequence |T]). But for finite n the 
calculation can be done explicitly by computer algebra. 
Lee et al. carried out such a calculation up to n = 9 
fl8j |. We succeeded in reaching n — 14 by the help of the 
following recursion. 

Clearly, we know from the expansion of a single expo- 
nential that 



1 



IN 



r(*i-*o) 



\m\\ 



(82) 



This serves as starting point of our recursion which relies 
on 

U p+1 {t) =e -^ ( -i)p+i(Vh2-Vn)* Up (Q 



E 

weB 



{(-it)INI[(_i)P+V z ]Hi£E±L 



■4., 



A A \ 



(83a) 



x {-it) ll - l] (T l ™ ] C^A r , 



,A m2 A mi . (83b) 



meB 



The comparison of the arising coefficients with those in 
(f8"Tj) leads to the recursion relation 

f_-n(p+i)|«i| 

° V p + i= E Hffllll (84) 



(u;,m)— v 



where the sum over (w, to) = u means that all splittings 
of the binary word v in two subwords w for the first part 
and to for the second part are considered. Given v with 
letters there are \ \v\ \ + 1 such splittings. 
The recursion (|84[) can be easily implemented in com- 
puter algebra programmes such as MAPLE. With about 
2 Gigabyte RAM the verification of the vanishing of the 
Cn with odd checksum \v\ for the optimum sequence (fTJ) 
up to the order n = 14 was feasible. Nevertheless, a 
general mathematical proof would be highly desirable. 



better than the concatenated sequence (CDD) proposed 
by Khodjasteh and Lidar [U, El]. [Ill But they found 
that qubits made from phosphorous impurities in silicon 
are better dynamically decoupled by the CDD sequence. 
They relate this result to the applicability of an expan- 
sion in time. Their model consists of a qubit coupled to a 
spin bath so that a direct applicability of results obtained 
for the spin-boson model is not possible. Yet the ques- 
tion is intriguing whether one can mimic the qualitative 
aspects of the spin bath by a bosonic bath. 

From the way the general single axis model is treated 
to derive the effect of the UDD sequence, see previous 
section, it is clear that the expansion in powers of t plays 
the crucial role. If such an expansion in time does not 
work, for instance because the resulting expansion is only 
asymptotically valid, there is no justification to use the 
UDD sequence. 

The analytically accessible spin-boson model allows us 
to investigate the question of the expansion in time in 
a concrete example. Inspecting Eqs. (|40b|43b|45d|) one 
realizes that the existence of the expansion of the signal 
s(t) depends on the existence of the expansions of Xn(t) 
and (f n (t)- In order that Xn(t) = 0(t n+1 ) the first n 
derivatives of Xn(t) must exist and vanish and the n + 
1st derivative must exist. From (|43b[) we see that the 
expansion of Xn(t) in powers of t is directly related to 
the expansion of y n (z) as in (|42[) in powers of z since 
z = cot. In Subsect. Ill CI and in Ref. [l6| we considered 
only the existence and the vanishing of the derivatives of 
y n (z). The existence of the integral over the frequencies 
is no issue as long as a completely hard cutoff at lu-£> is 
considered 



Joo(w) := 2aujQ(ujY) — u>) 



(85) 



for which no ultraviolet (UV) divergence can appear. 
Hence all derivatives with respect to time exist for Xn(t) 
and for (p n (t). This remains also true if the UV cutoff is 
exponential. 

But the physical systems might be such that the UV 
cutoff is soft because the spectral density displays power 
law behavior. We consider 



JJuj) :-- 



2au> 



1 + (UJ / WD )T 



(86) 



as generic form for this situation. Note that 7 = 00 
amounts up to the completely hard cutoff. The vanish- 
ing of the first n derivatives of y n (z) implies y n {z) = 
A[uit) n+1 plus higher terms. But in order to be able to 
conclude that Xn(t) = Ct ln+2 plus higher terms the in- 
tegral 



IV. INFLUENCE OF THE HIGH-ENERGY 
CUTOFF 



A 2 r°° 

C=— uj 2n J(uj)coth(Poj/2)duj (87) 
4 Jo 



Lee et al. [l8( observed that the optimized sequence 
((T|), henceforth abbreviated UDD, works very well in nu- 
merical simulations for GaAs quantum dots where it does 



must exist, i.e., converge. For J 7 (w) this strictly requires 
7 > 2n + 2. (88) 
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The equivalent consideration for the phase ip n (t) leads 
to a less strict condition. If y n (z) — A(ujt) n+1 plus higher 
terms one has ip n (t) = Dt n+1 which contributes the same 
order t 2n+2 as exp(— 2x n ) to 1 — s(t) because of the cosine 
in which it appears, see (|40b|) . The coefficient D is given 
by the integral 



D = — — I lu"- 1 J(w)du. 



Its existence requires only 



7 > n + 1 



(89) 



(90) 



for J 7 (w). Hence we conclude that the condition for the 
smallness of the deviations resulting from \n implies the 
condition for the smallness of the deviations resulting 
from (p n . For this reason, we will focus on the condition 
for the smallness of the exponential suppression by % n . 

For practical purposes, the existence or non-existence 
of certain derivatives is not the ultimate criterion. So 
below we compare the effect of various pulse sequences on 
the signal s(t). First, we look at the sequence UDD„(t) 
which is characterized by |T]). It leads via (|4"2")) for an 
even number of pulses n to 



y™°{z) = -2ie^ 2 sin(|) + 



^/2 



2 £(-1)^ sin (-cos 




Recall Eq. ([76]) stating the order y^ DD {z) = <D(z n+1 ). 
Second, we consider the concatenated sequence (CDD) 

[HE!]. The zeroth level CDD oW is the evolution with- 
out pulse. Higher levels are defined recursively by 

CDD l+1 (t) = CDD;(^) o CDD,(^) V/ odd (92a) 

CDD i+ i(t) = CDD^jonoCDDi^) VZ even, (92b) 

where o stands for the concatenation and n for a tt pulse. 
We obtain for the CDD sequence 



Hi 



CDD 



(z) = {-2i) l+1 e lz / 2 sin(2-'- 1 ^) \\ sin^-^z), 



k=l 



where I now stands for the level which is exponentially 
related to the number of pulses n w 2 l . From (|93|) it is 
easy to see that yf DD (z) = 0(z l+1 ) holds. 

Third, we consider the first suggestion 0, 0] , namely 
the periodic bang-bang (BB) control with n pulses and 

Sj=j/(n+l) (94) 

implying (for even n) 

y® B (z) = -2ie lz ' 2 cos (z/2) tan (z/(2n + 2)) . (95) 



From this equation one learns y BB (z) = O(z). 

Fourth, we consider the Carr-Purcell-Meiboom-Gill 
(CPMG) sequence @, H, 0] • This sequence results from 
the fc-fold iteration of a two-pulse cycle of length r = t/k . 
The pulses occur at t/4 and 3t/4. This cycle corresponds 
in fact to UDD2 (r) [l6| . We will come back later to iter- 
ations of UDD sequences. Here we state that CPMG is 
characterized by 



Si = (j - l/2)/n 
implying (for even n) 

y% FMG (z)=4ie iz / 2 sm 2 {z/(4n)y 



(96) 



sin(z/2) 



cos(z/(2ri)) 

From this equation it is clear that y^ PMG (z) — 0(z 3 ) 



(97) 



t- 10 




FIG. 1: Various pulse sequences (see main text) are compared 
for various values of the cutoff parameter 7 in J T (w), see (I86|l . 
All sequences comprise n — 10 pulses; this values is chosen for 
better comparison because level I = 4 of the CDD sequence 
has 10 pulses. The coupling value a in the spectral densities 
J-t(u>) is fixed to 1/4 and the temperature is zero. 

In Fig. [TJ the four sequences are compared for 10 ir 
pulses at a fixed value a = 1/4 of the coupling to the 
bath. The results for other values of a are very simi- 
lar. Furthermore, the temperature is fixed to T — be- 
cause the precise value of the temperature matters only 
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for small frequencies to — > while we focus here on high 
frequencies and the UV cutoff. 

In all six panels it is obvious that the bang-bang (BB) 
sequence does worst in accordance with the power law 
which is only linear. This inefficient suppression of deco- 
herence also implies that phase effects in the signal s(t) 
due to (p n in Eq. (|40b[) ) are seen most strongly leading to 
the bumps in Fig. [TJ We conclude that one should always 
try to use one of the other sequences. 

Comparing the CDD and the CMPG sequences the 
CPMG sequence is almost everywhere advantageous. 
Only for very low deviations 1 — s(t) the CDD does better 
because its curve is steeper reflecting a higher order in t: 

yCDD = 0^5) while y CPMG = ( f 3y 

Comparing the CDD and the UDD sequences the UDD 
sequence yields always lower deviations, except for very 
soft cutoffs (7 = 2) where both sequences behave equally. 
We conclude that we cannot explain the behavior found 
by Lee et al. [HI where the CDD sequence seemed to out- 
perform the UDD on the basis of the spin-boson model. 
Note that the slope of both sequences in Fig. Q] seems 
to be similar though this is difficult to tell from the 
depicted range of parameters. But the analytic results 
clearly states y% BB (z) = 0(z 5 ) while ^ DD (z) = 0{z n ) 
for the same number of pulses, namely n = 10. 

The interesting issue is a comparison of the CPMG 
and the UDD sequence. For very soft cutoffs, i.e., low 
values of 7, the CPMG sequence is slightly better. This 
was also observed in a model of classical Gaussian noise 
[3II ]. The UDD sequence, however, performs better for 
large values of 7. Indeed, this finding supports our an- 
alytical condition ([88]). As long as 7 ^ 2n the CPMG 
sequence with its relatively low order t 3 (in y n (uit)) does 
slightly better than the high-order UDD with i 11 . But for 
7 gl 2n the UDD outperforms the CPMG, especially for 
low deviations 1 — s(t) which matter most for quantum 
information processing. 

We substantiate the comparison between the UDD and 
the CPMG sequence further by Fig. [2] The results go 
into the same direction as before. As long as n ;$ 7/2 
the UDD does significantly better than the CPMG. For 
n w 7 the UDD does better than the CPMG at low values 
of 1 — s{t) ss 10~ 4 while the CPMG is advantageous at 
higher values 1 - s(t) w 10 _1 . For n > 7, the CPMG 
does slightly better than the UDD except for very small 
values of 1 — s(t). This constitutes a clear message for 
applications. 

One may wonder whether there is a way to combine 
the advantages of the UDD and of the CPMG sequence. 
Indeed, this is possible by resorting to hybrid solutions 
proposed already earlier [ijj ■ The UDD cycles with 
low values of n can be iterated. We denote such a se- 
quence by iUDD miC (i) where m stands for the number of 
pulses within one cycle and c for the number of cycles so 
that n = mc is the total number of pulses. This means 
we consider the concatenation 

iUDD ro>c (t) = (UDD m (t/c)) c . (98) 




FIG. 2: Comparison of the performance of the CPMG (dashed 
lines) and the UDD sequence (solid lines) for various number 
of pulses n (legend holds for dashed and solid lines) and a 
spectral density with 7 = 8 corresponding to a cutoff of in- 
termediate hardness. Other values are fixed: a — 1/4, T — 

A quantitative comparison for iterated iUDD se- 
quences is given in Fig. [3] for a total of 12 ir pulses. Note 
that iUDD 2 ,e is equivalent to CPMG while iUDDi 2 ,i is 
equivalent to the UDD sequence. The guideline here is 
the corollary 

7 > 2m + 2 (99) 

of |88|) where m is the number of pules within one cycle. 
It results from the observation that Vm^i 2 ) is of order 
z m+i i n d e p en d en t from the number of cycles. 

If the condition (|M)l is not valid the use of any sequence 
of higher order does not pay. This is clearly seen in the 
uppermost panel for 7 = 4 (very soft cutoff) in Fig.[3j All 
curves are almost on top of each other. The CPMG, i.e., 
iUDD2,6, is slightly better than the other pulse sequences. 

In the middle panel for 7 = 8 (intermediate cutoff) 
in Fig. [3J the situation has changed. For low deviations 
1 — s(t) the use of the iUDDa^ or the iUDD^ sequence 
pays while the implementation of a larger value of m 
hardly pays. 

In the lowermost panel for 7 — 00 (hard cutoff) in Fig. 
[3] the implementation of higher order sequences is always 
useful for low values of 1 — s(t) as was to be expected. 

Fig. [3] illustrates that one can gain considerably in co- 
herence without implementing the fully optimized pulse 
sequence (|T|). Already the implementation of periodic cy- 
cles with a moderate number of pulses can be very help- 
ful. In practice, this strategy is generally much easier to 
realize since not so many special instants in time need to 
be fine-tuned. 

Another remark for experimental realizations is in or- 
der. If the pulses are not ideally tailored then the advan- 
tageous of dynamical decoupling will be thwarted by ac- 
cumulated pulse errors. So in practice one always will be 
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FIG. 3: Comparison of the performance of various iterated 
iUDD mjC (t) sequences for a total number of n — 12 pulses and 
various cutoff exponents 7. Other values are a — 1/4, T = 0. 

faced with the need to find the optimum tradeoff. Note, 
however, that this fact makes it particularly interesting 
to reach an optimum suppression of decoherence with a 
small number of pulses. 

V. CONCLUSIONS 

We investigated the suppression of decoherence by se- 
quences of ideal, instantaneous 7r pulses. The model un- 



der study is a spin-boson model valid for pure dephasing, 
i.e., for a finite T 2 but an infinite T±. But also the most 
general model for phase decoherence (single-axis deco- 
herence) is considered. 

First, we have provided the detailed derivation of the 
equations which were used in our previous Letter in Ref. 
Ilq . In particular it was rigorously shown that the se- 
quence (JTJ) (UDD) makes the first n derivatives vanish. 
Furthermore, it was shown that the results transfer also 
to the classical case of Gaussian fluctuations. 

Second, it was shown that the UDD sequence is advan- 
tageous for any initial state. This important finding was 
achieved by analyzing the corresponding time evolution 
operator. 

Third, we considered the most general model for phase 
decoherence and extended the analytical results of Lee 
et al. to the 14th order in the time. This was achieved 
on the basis of an efficient recursion scheme suitable for 
implementation in a computer algebra programme. 

Fourth, we investigated the influence of the high- 
energy cutoff in the framework of the single-axis spin- 
boson model. We compare various pulse sequences 
which are currently under debate, namely the periodic 
bang-bang sequence (BB), the concatenated dynami- 
cal decoupling (CDD), the well-established Carr-Purcell- 
Meiboom-Gill sequence (CPMG) and the general itera- 
tion of UDD cycles (iUDD). 

The most important observation is that decoherence 
due to baths with very soft cutoffs are much more dif- 
ficult to suppress than decoherence due to baths with 
hard cutoffs. For soft cutoffs, the simpler sequences 
(CPMG= iUDD 2 , c or iUDD mjC with low values of m) 
are completely sufficient. Higher order sequences do not 
pay. We established a rule of thumb when the implemen- 
tation of a more intricate sequence is appropriate. The 
number of pulses m in one cycle should not exceed 7/2 
where 7 is the exponent of the high-energy power law of 
the decohering spectral density </ 7 (w), see Eq. (I5FJ)) . 

By the above results, we have elucicated the possibili- 
ties of dynamical decoupling. Mathematically, important 
derivations are provided. Practically, important guide- 
lines are established under which conditions which se- 
quences are most appropriate. 



Acknowledgments 

I like to thank J. Cardy, D. Lidar, S. Pasini, P. Kar- 
bach, T. Fischer, J. Stolze, D. Suter, L. Viola, and W. 
Witzel for helpful discussions. 



12 



[1] E. L. Hahn, Phys. Rev. 80, 580 (1950). 
[2] H. Y. Carr and E. M. Purcell, Phys. Rev. 94, 630 (1954). 
[3] S. Meiboom and D. Gill, Rev. Sci. Inst. 29, 688 (1958). 
[4] U. Haeberlen, High Resolution NMR in Solids: Selective 

Averaging (Academic Press, New York, 1976). 
[5] P. Zoller, T. Beth, D. Binosi, R. Blatt, H. Briegel, 

D. Bruss, T. Calarco, J. I. Cirac, D. Deutsch, J. Eisert, 

et al., Eur. Phys. J. D 36, 203 (2005). 
[6] L. Viola and S. Lloyd, Phys. Rev. A 58, 2733 (1998). 
[7] M. Ban, J. Mod. Opt. 45, 2315 (1998). 
[8] L. Viola, E. Knill, and S. Lloyd, Phys. Rev. Lett. 82, 

2417 (1999). 
[9] L. Viola, Phys. Rev. A 66, 012307 (2002). 
[10] L. F. Santos and L. Viola, arXiv:0804.0890. 
[11] W. Yao, R. B. Liu, and L. J. Sham, Phys. Rev. Lett. 98, 

077602 (2007). 

[12] W. M. Witzel and S. Das Sarma, Phys. Rev. Lett. 98, 
077601 (2007). 

[13] K. Khodjasteh and D. A. Lidar, Phys. Rev. Lett. 95, 
180501 (2005). 

[14] K. Khodjasteh and D. A. Lidar, Phys. Rev. A 75, 062310 
(2007). 

[15] W. M. Witzel and S. Das Sarma, Phys. Rev. B 76, 241303 
(2007). 

[16] G. S. Uhrig, Phys. Rev. Lett. 98, 100504 (2007). 

[17] D. Dhar, L. K. Grover, and S. M. Roy, Phys. Rev. Lett. 

96, 100405 (2006). 
[18] B. Lee, W. M. Witzel, and S. Das Sarma, Phys. Rev. 

Lett. 100, 160505 (2008). 
[19] H. G. Krojanski and D. Suter, Phys. Rev. Lett. 97, 

150503 (2006). 

[20] L. M. K. Vandersypen and I. L. Chuang, Rev. Mod. Phys. 

76, 1037 (2004). 
[21] E. Fraval, M. J. Sellars, and J. J. Longdell, Phys. Rev. 

Lett. 95, 030506 (2005). 
[22] J. R. Petta, A. C. Johnson, J. M. Taylor, E. A. Laird, 



A. Yacoby, M. D. Lukin, C. M. Markus, M. P. Hanson, 
and A. C. Gossard, Science 309, 2180 (2005). 

[23] J. J. L. Morton, A. M. Tyryshkin, A. Ardavan, S. C. Ben- 
jamin, K. Porfyrakis, S. A. Lyon, and G. A. D. Briggs, 
Nature Phys. 2, 40 (2006). 

[24] A. Greilich, R. Oulton, E. A. Zkukov, I. A. Yugova, 
D. R. Yakovlev, M. Bayer, A. Shabaev, A. L. Efros, I. A. 
Merkulov, V. Stavarache, et al., Phys. Rev. Lett. 96, 
227401 (2006). 

[25] A. Greilich, , D. R. Yakovlev, A. Shabaev, A. L. Efros, 
I. A. Yugova, R. Oulton, V. Stavarache, D. Reuter, 
A. Wieck, et al., Science 313, 341 (2006). 

[26] S. Pasini, T. Fischer, P. Karbach, and G. S. Uhrig, Phys. 
Rev. A 77, 032315 (2008; S. Pasini and G. S. Uhrig, J. 
Phys. A: Math. Theor. 41, 312005 (2008 ). 

[27] A. J. Leggett, S. Chakravarty, A. T. Dorsey, M. P. A. 
Fisher, A. Garg, and W. Zwerger, Rev. Mod. Phys. 59, 
1 (1987). 

[28] U. Weiss, Quantum Dissipative Systems (World Scien- 
tific, Singapore, 1999), 2nd ed. 

[29] D. Li, A. E. Dementyev, Y. Dong, R. G. Ramos, and 
S. E. Barrett, Phys. Rev. Lett. 98, 190401 (2007). 

[30] D. Li, Y. Dong, R. G. Ramos, J. D. Murray, K. MacLean, 
A. E. Dementyev, and S. E. Barrett, Phys. Rev. B 77, 
214306 (2008). 

[31] L. Cywihski, R. M. Lutchyn, C. P. Nave, and S. Das 
Sarma, Phys. Rev. B 77, 174509 (2008). 

[32] Note that we have changed the definition of f(t) by a 
factor of 2 relative to our previous work [l6l ] to keep the 
notation concise. 

[33] It should be mentioned that the CDD was originally pro- 
posed for the general decoherence model where all Pauli 
matrices of the qubit are coupled to the bath. In the 
present article, we consider the special case of CDD for 
single axis models for phase decoherence. 



